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ABSTRACT 

A non-similar boundary layer theory for air blowing over a water layer on a flat plate 
is formulated and studied as a two-fluid problem in which the position of the interface is 
unknown. The problem is considered at large Reynolds number (based on x), away from 
the leading edge. We derive a simple non-similar analytic solution of the problem for which 
the interface height is proportional to x 1 / 4 and the water and air flow satisfy the Blasius 
boundary layer equations, with a linear profile in the water and a Blasius profile in the air. 
Numerical studies of the initial value problem suggests that this asymptotic, non-similar 
air-water boundary layer solution is a global attractor for all initial conditions. 
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1 Introduction 


The effects of water layers driven over solid surfaces by wind are of interest in the performance 
of aircraft in rain, for the de-icing of airplane wings and surely in many other applications. 
Since such problems are intrinsically of a boundary layer type and since the more interesting 
phenomenon which might arise, like the formation of waves, film rupture and the like are 
probably best framed in terms of stability, it is necessary to derive the analytic forms that 
such flows will take on when instability is neglected. This derivation is carried out here. 

Previous works related to the present one are by Yih (1990), who modeled the de-icing 
problem, and by Wang (1992), who considered the development of boundary layers in the 
shearing flow of one fluid over another. Both works are flawed by assuming rather than 
finding the shape of the interface. In the case treated by Wang, the interface is flat but the 
jump of the normal stress is not zero. Yih assumed that flow in the air is a Blasius flow and 
the flow in the water is a simple shear. This is a correct form for the boundary layer, as we 
shall show, but he neglects the variation of film thickness with x and so cannot enforce the 
kinematic condition at the interface or find its shape. 

In this work, we formulate a non-similar boundary layer theory retaining all terms which 
decay faster than l/£, £ = x(U/2iy 2 x) 1 ' 2 , where v 2 is the kinematic viscosity of air. All of 
the interface conditions are enforced in this asymptotic regime. An effect of the small but 
non-zero vertical velocity component at the interface is to force the interface to grow like 
x x ! A when the boundary layer in the air grows like x ^ 2 . The interface looks thin on the 
scale of the boundary layer. Asymptotically, at large x, the water and air satisfy the Blasius 
boundary layer equations with a linear profile in the water and the flat plate profile in the 
air. This non-similar (or coupled self-similar) solution appears to be a global attractor for 
all initial conditions. 


2 Governing Equations 

A water film of height y = h(x) is flowing on a flat plate driven by shear stresses emanating 
from an air stream with streaming velocity U. We seek the nature of the flow under the 
circumstances which give rise to Blasius boundary layers in the flow of one fluid over a flat 
plate. This motivates the introduction of the same scales that are used in the classical case, 
giving 
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where u = /' is the x component of velocity, v = p}' — f — £ §£ is the y component if velocity, 
t is the time, p is the pressure and / is the stream function. The parameters p 2 and v-i are 
for air; subscript 1 is for water. 

These scales are introduced into the Navier-Stokes equations which are written below in 
the new variables without approximation. The continuity equation becomes 
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I 11 the interface equations written below, [.] designates the jump (.)i — (.) 2 - At the interface 
p = p*(y = h(x)), the velocity is continuous, 
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The shear stress is also continuous across p = p*, 
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The jump in the normal stress is balanced by interfacial tension, 
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where G = (g/U 2 ) {2v^ x/U) 1 ^ 2 , g is the acceleration due to gravity, S = T /(V2P2U), and T is 
the coefficient of surface tension. Density and viscosity are normalized by the corresponding 
properties of the air 
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The kinematic equation of the free surface is written in the boundary layer coordinates 
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3 Asymptotic solution for large £ 

Here we exhibit a special solution of our coupled air-water system which will be shown later 
to be a global attractor of non-similar boundary layer solutions. We start by making several 
assumptions. First we assume that for large x (large £) the horizontal velocity component 
in the water is linear in y, and there is a similarity solution in the air with /( 0) = /'( 0) = 0, 
and /"( 0) is constant. In the air we write 
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Secondly we assume that the interface position 77* — ► 0 as £ — » 00. This assumption is 
equivalent to assuming that as £ — > 00 the liquid layer is a vanishingly small fraction of the 
boundary layer in the air. 

We are going to show that the continuity of the shear stress across the interface 77* implies 
that c = Ax* 1 / 2 , with A a constant, so that u = At) can be expressed in terms of 77 alone. 
The interface is on y = h(x),r) = 77* = k-£jj. The shear stress is continuous 
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Now for large x, we have assumed that 77* = h(x)/x 1 / 2 — » 0, so that f"(r)*) — > /"( 0) which 
is constant. Then 
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Using (13), we may write 
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Using (18), we find that 
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Now 

f(v*) = u(h(x),x). 

For large x,ri* -* 0, so that /'(»/*) — ► /'( 0) = 0. Actually we know that u tends to zero 
because the mass flux is constant but the flow area (18) keeps increasing. We may also show 
that v(h(x),x) tends to zero for large x. We can do this in different ways. Suppose we 
multiply the kinematic condition (11) by £u and let £ — * oo. 
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From (19) we know that, since £ = kx 1 ^ 2 , as £ — * oo, 
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Since u — ♦ 0 as £ — * oo, so also v — ► 0. So as £ — * oo, both u(h(x), x) — * 0 and v(h(x), x) —* 0; 
that is — > 0 and /(»/*) — > 0. The interface conditions then become the same as the 

boundary conditions for a Blasius boundary layer in the absence of the water layer but 
displaced, and so the velocity in the air is given by the Blasius solution for the boundary 
layer over a flat plate. The solution in the water f = Arj also satisfies the Blasius equation 
//" + V2f"' = 0, 0 < t} < ri* —> 0 and implies that the jump conditions [/] = [/'] = 0 may 
be replaced with the Blasius boundary conditions /( 0) = f{ 0) = 0 applied upon a solid 
surface. 

This shows that our original assumptions are self-consistent. The reader may^ verify 
that the solution (13) with c(x) given by (16) is self-similar (satisfies (24) without the the 
bracketed term). The main good luck that we have is that 77 * — »■ 0 for large x even though^ 
h(x) — ► 00 
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4 Non- similar boundary layers 

Since the asymptotic solution for large £ is similar, but not in the usual sense, we are 
motivated to see if this solution can be embedded in a large class of non-similar solutions for 
large £ (large Reynolds number Uxj 2v 2 ) and to investigate the possibility that the asymptotic 
solution enjoys a special status as a global attractor for all solutions of this non-similar family. 

Self similar boundary layers depend only on 7 / and not on £. The £ derivatives of u, v and 
Tf* may be small, of the order of £~ n ,n > 1, when £ is large. We retain those terms in each 
equation which are 0(1), dropping all terms of 0(£ -n ), For this computation we assumed 
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that the first derivative of u with respect to £ scale with £ -1 , the second derivative with £ 2 
and = 0(£~ 3/2 ). These scalings can be verified a posteriori. We find that 


ff" + „/» _ = o, 


(24) 


and at tj = 77 *: 

1 = 0, [/1 = 0. 1/1 = 0, [Rf] = 0 , M = = -’'■)• < 25 ) 

Equation (24) can be found in Schlichting (1987). We may remark that the contributions 
of the viscous terms in the normal stress balance are £ -2 times the terms retained and the 
surface tension terms are £ -3 times those retained. 
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5 Solutions of initial value problems for the non-similar, 
two-fluid boundary layer equation 

5.1 Evolution of a profile which is initially parabolic in segments. 

Since the non-similar solutions depend on £, we are obliged to consider the evolution of flows 
prescribed at some initial position £ = Fortunately our simulations of those initial value 
problems indicate a rapid decay from initial values to the asymptotic solution described in 
§3. This situation is not unfamiliar. In the classical theory of boundary layers, the Blasius 
solution, the similarity solution, is a global attractor for all initial value problems which are 
not similar (see Serrin (1967)). In our case the attractor cannot be self-similar, but the non- 
similar solution of §3 arises asymptotically, for large £ from all initial conditions explored. 
The first initial condition we used to solve (24) and (25) is plotted below. 



(a) (b) 

Figure 1. Parabolic initial condition assumed for the horizontal velocity component /' with 
rj* — 0.15 and £ 0 = 50. (a) Initial condition, (b) Detailed plot of initial condition at the 
interface. 

For this initial condition r/* = 0.15 and £ 0 = 50. These values were chosen to be representative 
of conditions discussed by Hastings and Manuel (1985), in which they describe the results 
of their wind tunnel measurements of a wing in simulated rain. In the above profile, for 
0 < T] < T}*, f is parabolic with /'" being positive. For rj* < rj < 3.5, /' is parabolic with f" 
being negative, and for r] > 3.5, f = 1.0. 

We solved (24) and (25) subject to the aformentioned initial conditions for £ > £ 0 using 
a finite difference scheme found in Schlichting (1987). This scheme is iterative using second- 
order differencing for the rj derivatives, and a first-order, implicit Euler differencing for the £ 
derivatives. We assumed an interface rj* = a£ n , where a = rj* (£„) / an< ^ found the n(£) 
given in figure 2. 
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Figure 2. Plot of n for the non-similar boundary layer with initial conditions from figure 1. 

The exponent function n(£) decays rapidly to -1/2, corresponding to a limiting power 
law increase h(x ) oc x 1 / 4 of the interface height as given by the asymptotic solution shown 
in figure 1. For £ not much larger than £ 0 the horizontal velocity component f becomes 
linear in rj, f = c(()rj, in the water layer and remains so for all greater £. The shear stress 
in the water, j£c, is directly related to the shear stress in the air at the interface through 
the tangential stress condition found in (25). From our calculations we find that the shear 
stress in the air at the interface, f", asymptotically tends to the shear stress in the Blasius 
boundary layer of a single fluid over a flat plate at the plate surface, as shown in figure 3. 




Figure 3 (a) Detailed plot of boundary layer solution in the water layer at £ = 50.0675 for 
the initial conditions of figure 1. (b) Shear stress, f", in the air at the interface, jj* = 


for £ > &>. 


Notice that the exponent function n(£) tends to its asymptotic value extremely rapid, while 
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the shear stress at the interface in the air f" has a much slower decay. The reason for this 
difference is n( £) becomes nearly -1/2 as soon as the horizontal velocity component in the 
water becomes linear, which occurs for ( not very much greater than ( 0 . Figure 3 shows that 
the flow has not yet reached its asymptotic form, since the shear stress at the interface in 
the air is appreciably different from its asymptotic value. The flow then slowly tends to its 
asymptotic state. The flow is non-similar for finite £ because the boundary layer height and 
the interface do not scale the same. Asymptotically one may say the flow is similar because 
rj* goes to zero and the flow in the air is the self-similar Blasius profile. 

5.2 Evolution of another profile which is initially parabolic in 
segments. 

In the next example we used the same initial /'(tj) for t] > rj* as in figure 1(a), but the profile 
shown in figure 1(b) was replaced with the one shown in figure 4. 



Figure 4. Initial profile in the water at £ = £ 0 . The profile in the air is shown in figure 1(a). 

In figure 5 we show the evolution of the exponent n(£) for £ > £ 0 for the aformentioned 
initial condition. Obviously n(£) — ► —1/2 for large £ and r)* = a£ -1 / 2 for large £. The 
solution given in §3 is attained asymptotically. 
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5.3 Evolution of a profile which is initially in similarity form but 
with 77 * (£ 0 ) = 0.15 > 0. 

The initial condition shown in figure 6 is generated from (24) and (25) when the ( derivatives 
are set to zero, so that 

SI" + »!"' = 0, (26) 

and at r) = t)*: 

|^ = 0, 1/1 = 0, [/] = 0, [*/"] = 0, \p) = [pW. (27) 

This system is solved using a iterative finite difference scheme using second order differencing 
and gives rise to the profiles shown in figure 6. 



(a) (b) 

Figure 6. Solution of (26) and (27). (a) 77 > 0.15. (b) t ] < 0.15. The Blasius solution 
of (24) for /(0.15) = /'(0.15) = 0 is also shown in (b). 



9 


The profiles in figure 6 were used for initial values in the non-similar equations (24) and 
(25) and they evolved to the coupled self-similar solutions described in §3 in which — ► 0. 

The difference between the solution in the air and the classical Blasius flat plate solution 
also tends to zero. A graph of the exponent function for this example is shown figure 7. 



Figure 7. Exponent function n(£) for rj* = for the initial conditions of figure 6. 

In these examples and in all the others, which we tried but are not shown, the asymptotic 
solution given in §3 is ultimately attained. 


6 Behavior of £ derivatives of the non-similar solution 

In deriving the system which lead to the above solutions, we assumed that first and second 
derivatives of u and r\* with respect to £ were inversely proportional to £~ n ,n > 1. ^From 
our solution we find that %- scales like £~ 3 ' 2 and scales like ( ~ 5 / 2 , in agreement with 
our scalings. Since £ scales with a: 1 / 2 , the interface position scales with x 1 / 4 . 

In order to examine the £ derivatives of u, we define two functions, 

F,(l,()= (28) 

( 29 ) 

Depending upon the initial conditions chosen at £ 0 , Fi(rj,£ 0 ) and F 2 (r/,£ 0 ) may be extremely 
large. However, for larger values of ^ the horizontal velocity component in the water becomes 
linear in 77, n(£) tends to -1/2 and F\(r},() < 1 , < 1 (see figure 8). 
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Figure 8. Plot of Fi(t), £) and F 2 (tj() versus rj for £ = 50.01 for the flow plotted in figure 6. 

W and \F 2 \ decrease, as shown in figure 9, although the decrease in |Fj| is very small. 
Since as ( — ♦ oo, »?* — » 0 and the flow in the air goes to the Blasius boundary, both F\ — * 0 
and F 2 — » 0 as £ — ► oo. 




Figure 9. Plot of and F 2 (i]Q versus rj for £ = 129.1 for the flow plotted in figure 6. 

Then since F\ and F 2 are both 0(1), the assumptions that were made in deriving (24) and 
(25) are shown to hold. 
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7 Concluding remarks 

The analysis given in this paper can be extended to two-fluid boundary layer problems with 
other free streams, say U — UoqX 71 . It is also probable that the solution given in §3 and the 
other limiting solutions to which we have just alluded are unique large x limits of steady 
coupled air-water solutions of the Navier-Stokes equations with different initial profiles at 
x = x 0 . 
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